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ABSTRACT 

Hot super-Earths likely possess minimal atmospheres established through va- 
por saturation equilibrium with the ground. We solve the hydrodynamics of 
these tenuous atmospheres at the surface of Corot-7b, Kepler 10b and 55 Cnc-e, 
including idealized treatments of magnetic drag and ohmic dissipation. We find 
that atmospheric pressures remain close to their local saturation values in all 
cases. Despite the emergence of strongly supersonic winds which carry sublimat- 
ing mass away from the substellar point, the atmospheres do not extend much 
beyond the day-night terminators. Ground temperatures, which determine the 
planetary thermal (infrared) signature, are largely unaffected by exchanges with 
the atmosphere and thus follow the effective irradiation pattern. Atmospheric 
temperatures, however, which control cloud condensation and thus albedo prop- 
erties, can deviate substantially from the irradiation pattern. Magnetic drag and 
ohmic dissipation can also strongly impact the atmospheric behavior, depend- 
ing on atmospheric composition and the planetary magnetic field strength. We 
conclude that hot super-Earths could exhibit interesting signatures in reflection 
(and possibly in emission) which would trace a combination of their ground, 
atmospheric and magnetic properties. 



1. Introduction 

Hot super-Earths are an emerging class of exoplanets. CoRot 7b, Kepler 10b and 55 
Cnc-e are currently the only well-characterized, transiting members of this new class (Leger 
et al. 2009; Batalha et al. 2011; Fischer et al. 2008; Winn et al. 2011; Demory et al. 2011) 
but observational trends suggest that many more hot super-Earths will be discovered and 
characterized in the future (e.g., Howard et al. 2011; Borucki et al. 2011). 

Structural and evolutionary considerations suggest that hot super-Earths are rocky plan- 
ets which have lost their original volatile atmospheres, and possibly even some of their rocky 
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material, via atmospheric erosion (Valencia et al. 2010). Indeed, a hot rocky super-Earth 
will retain a minimal atmosphere that is continuously replenished via vapor saturation equi- 
librium with the ground (Schaefer & Fegley 2009; Leger et al. 2011) and is thus perpetually 
subject to erosion. Hot super-Earths are also expected to be tidally-locked to their parent 
stars. Under these conditions, large temperature and surface pressure differences will exist 
between the day and the night sides of these planets. By contrast with the majority of 
known planetary atmospheres, which have a well distributed atmospheric mass around the 
planet, an unusual circulation regime will thus develop at the surface of hot super-Earths, 
with powerful winds dynamically redistributing atmospheric mass via exchanges with the 
ground. 

Understanding the atmospheric behavior of hot super-Earths is important for the in- 
terpretation of direct observational data and because it informs evolutionary considerations 
such as atmospheric erosion. In this Letter, we study the unusual hydrodynamics of the at- 
mosphere of hot super-Earths, using a simple model inspired by related work on lo (IngersoU 
et al. 1985). In §2, we describe our model hypotheses, equations and method of solution. 
Our main results are presented in §3. We conclude by discussing some potential implications 
of our work in 54. 



2. Model 
2.1. Hypotheses 

We consider hot super-Earths that are tidally-locked to their parent stars and thus 
possess permanent day and night sides. We treat the atmosphere as a continuous fluid that 
is hydrostatically bound to the planet. The fluid approximation is justified, even if only 
marginally so in the most tenuous regions described by our solutions, where one runs out 
of atmosphere near the planetary day-night terminator. For simplicity, we ignore rotation 
and Coriolis effects in our analysis, so that point symmetry around the substellar point can 
be assumed. Comparing the rotation timescale (1/^^) for the three super-Earths of interest 
with the advection time across a planetary radius (r/V) we find Rossby numbers ~ 2 for the 
typical velocities ~ km/s obtained in our solutions. It is thus reasonable to neglect rotation 
as a first approximation but improved models may need to account for it|l] 



^Rotation could break the substellar-point symmetry of our solutions and, under conditions of permanent 
hemispheric forcing, could lead to the formation of superrotating equatorial winds (Showman & Polvani 
2011). 
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The atmospheric composition that results from vapor saturation equihbrium with the 
ground on a hot super-Earth is a priori unknown. The analysis of Schaefer & Fegley (2009) 
for such sublimating atmospheres suggests that monatomic sodium is the dominant con- 
stituent on a hot super-Earth with a bulk silicate earth composition, prior to any fractional 
loss. According to their analysis, this is also the case with the largest overall atmospheric 
mass (surface pressures). Therefore, without loss of generality, we adopt a pure sodium 
composition in our models, so as to maximize atmospheric effects in the solutions, and we 
comment on the consequences of adopting other compositions. Since the atmospheres of 
interest are tenuous, typically 100 mbar or much less, we assume that the absorption optical 
depths in the thermal and the visible are negligibly small (^ 1). Non-radiative consid- 
erations, in particular dynamics and exchanges with the ground, will then determine the 
thermodynamical state of the atmosphere. 

We postulate, and verify a posteriori, that latent heat exchanges with the atmosphere 
have at most a small effect on the surface energy budget, which is dominated by radiative 
fluxes. We therefore assume that the ground temperature, T,, can be obtained by simple 
radiative balance for the permanent dayside: T, = (Tsub — Tas) cos^/^6' + Tag for 6* < 85 ° , 
where 6 is the angle away from the substellar point, Tg^b is the substellar temperature, and 
Tas is the antistellar temperature. For 85 ° < 6' < 110 ° , we assume that the temperature 
decreases linearly with 9 to Tas = 50 K, a value which may be reasonable on the basis of 
geothermal heating of the nightside (Leger et al. 2011), with little consequences on our 
results. The substellar surface temperature is evaluated simply as Tsub = T^\/ R* / D , where 
and i?^ are the stellar effective temperature and radius, D is the planet-star orbital 
separation and a negligible albedo is assumed. 

2.2. Equations 

Our formalism is directly inspired from that used by IngersoU et al. (1985) for the 
study of a frost atmosphere on lo. Indeed, we validated our implementation by reproducing 
the results for lo first, before moving on to super-Earths. Like IngersoU et al. (1985), we 
consider that strong vertical exchanges in the thin, turbulent atmosphere justify the use of 
vertically-integrated equations for the conservation of atmospheric mass, momentum and 
energy. Assuming that the velocity and the entropy per unit mass (for a dry adiabat) are 
constant with height in the atmosphere]^ the governing equations for the pressure, tempera- 



^Different assumptions for these vertical profiles would lead to qualitatively similar results, as shown 
explicitly by IngersoU et al. (1985) with a dry vs. wet adiabat comparison. We note, however, that these 
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ture and velocity at the base of the atmosphere can be written (see IngersoU et al. 1985 for 
details) : 

^ {VP sm9) = mE, (1) 



rg sin 9 dO 



^ [{V' + ^C,T) P sine] = —^JC,TP + r, (2) 



rg sin 9 d9 rg tan 9 



1 

rgi sin 9 d9 



V' 



2 



2 + CpT]VP sin9 



= (3) 



where r is the planetary radius and g is the planet's surface gravity. The three unknowns 
of the atmospheric flow are the velocity, V, the pressure, P, and the temperature, T, at the 
base of the atmosphere. The thermodynamic parameter (5 = TZ/{TZ + Cp), where TZ = k^/m 
is the gas constant and Cp is the specific heat at constant pressure. For sodium, with a mean 
molecular weight /i = 23 and a mass per atom m = iimu, we adopt the monatomic values 
Cp = 904 J K-i kg-i and f3 = 0.286. 

E, r, and Q are the rates per unit area of molecules, momentum and energy trans- 
ferred to the atmosphere by exchange with the ground. These fiuxes, which are mediated 
by a boundary layer of negligible vertical extent at the base of the atmosphere, are governed 
by sublimation/condensation processes, flow advcction, and turbulent exchanges by eddies. 
Our implementation of these surface-atmosphere boundary layer exchanges follows closely 
that adopted by IngersoU et al. (1985). The net flux of sublimating/condensating molecules 
is proportional to the difference between the vapor pressure P„ of the surface and the at- 
mospheric pressure P, so that E = {Ps - P)/{vsV2tt), where Vs = {kTs/my^ is the 
molecular thermal speed at the surface. We approximate the vapor pressure curve of sodium 
by fltting the sodium curve for the bulk silicate atmosphere model shown in Figure 1 of 
Schaefer & Fegley (2009) with the formula Ps{Ts) = Aexp{-B/Ts). We adopt A = 10^-^ Pa 
and B = 38, 000 K. 

The fluxes Q and r depend linearly on the quantities being transported, with the same 
transfer coefficients Wg and Wa- Energy gains from the surface arc proportional to the en- 
thalpy per unit of mass at the surface, CpTs, while losses scale as the sum of the kinetic 
energy and enthalpy per unit mass, + CpT). For momentum, gains from the surface 

are zero, while losses are proportional to V. This leads to the formulation 



assumptions could become invalid if sufficient stellar energy is deposited in the atmosphere for a stably 
stratified thermal profile to be established. 
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T = -PsWaV, (4) 

Q = PsWsCpT, - PsMVV'^ + CpT), (5) 

where Ps = mPs/ {kbTg) is the density in the boundary layer. Two different contributions 
are included in the definition of the transfer coefficients Wa and Wd- First, the contribution 
from the mean flow normal to the surface is proportional to E and is represented by the 
velocity — mE/pg. Second, the contribution from turbulent eddies is obtained from 
turbulent boundary layer results in the case E — We adopt the same expressions for 
the transfer coefficients Wa and as IngersoU et al. (1985). Finally, where needed, we 
use Sutherland's ideal gas formula to evaluate the dynamic viscosity of the atmospheric gas: 
7] = r]o{Ts/Tof'\To + C)/(T, + C), with r/o = 1.8 x IQ-^ kg m-^ s-\ To = 291 K and 
C = 120 K. We have verified that our results do not strongly depend on details of the 
viscosity law adopted. 

In sufficiently ionized regions, magnetic drag and ohmic dissipation can influence at- 
mospheric flows on hot super-Earths, if these planets possess strong enough magnetic fields 
(Gaidos et al. 2010; Tachinami et al. 2011; Driscoll & Olson 2011). Here, we build on the 
recent understanding gained in the context of hot giant exoplanet atmospheres (Batygin & 
Stevenson 2010; Perna et al. 2010a,b; Menou 2011). Assuming induced electric currents 
that are confined to the thin atmosphere of hot super-Earths, magnetic drag and ohmic 
dissipation may be included in the simplest possible way as additional, vertically-integrated 
contributions to the fiuxes of momentum and energy received by the atmosphere, through 
linear terms proportional to the inverse of a representative magnetic drag time, Tdrag- There- 
fore, in models including magnetic effects, we add contributions to r and Q which act to 
brake the winds and heat up the atmosphere according to: 



PV _ PV^ 

"^mag — ) Qmag — '^^^ • ("J 

fl'-' drag 9-'- drag 

Adopting a magnetic field strength i? = 1 G, detailed estimates of the resistivity {rj) 
based on the formalism described in Menou (2011) yield drag times Tdrag rj/B'^ that 
vary greatly over the planet for a sodium atmosphere. Drag times can be much shorter 
than a typical advection time, r /V ^ in the vicinity of the hot substellar point, but they can 
become negligibly large in poorly-ionized regions where T < 1000 K. Changing the assumed 
composition to a less easily ionized atom, e.g. from Na to Mg (another possibly abundant 
constituent; Schaefer & Fegley 2009), can also lead to dramatically reduced magnetic effects. 
Therefore, rather than exploring magnetic effects in detail here, we illustrate their conse- 
quences below by presenting two idealized models for Kepler 10-b with spatially uniform 
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values of T^rag = 10^ and 5 x 10^ s, which have been chosen to match the typical advection 
time in our solutions with velocities ~ km/s. 



2.3. Method of Solution 

We solve our set of differential and auxiliary equations for the unknowns P, T, V with 
a relaxation method (specifically, the algorithm described in Hameury et al. 1998). Special 
care must be taken when integrating these equations because the solution goes through a 
sonic point at an unspecified angle, Os- To find Og, we define ^ as a fourth variable on 
a numerical domain describing the subsonic region only and we apply physically sensible 
boundary conditions on this domain. 

At the substellar point, V{9 = 0) = is imposed for symmetry and T{d = 0) = Tg^b. 
The transonic solution must have a well-behaved value of dV/dO at the sonic point, which 
provides the two additional boundary conditions needed on the subsonic domain. Noting 
that the following first-order differential equation is satisfied by the velocity, 

^ V'^\dV /I 1\ ,,2^ rgVr ^ C„TPV 

requiring a smooth sonic transition leads to the two boundary conditions at ^ = 

n T^fQ ^ l\mE\ 

Our initial guesses for the relaxation method are based on the analytical solutions of 
IngersoU (1989), adapted to our super-Earth problem. Once 9s and the solution on the 
subsonic domain are known, using a small extrapolation across the sonic point, the integra- 
tion can proceed in the supersonic domain using a simple Euler method. The integration is 
stopped when the pressure is so low that the fluid approximation ceases to be valid. 
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3. Results 

We solve the hydrodynamics of sodium atmospheres at the surface of CoRot-7b, Kepler 
10b and 55 Cnc-e (see Table 1 for the specific parameters adoptedj^l- Figures 1-4 show the 
pressure, temperature, velocity and Mach number obtained for the base of the atmosphere, 
as a function of the angular distance from the substellar point. In each figure, large dots 
indicate the location of the sonic points. Solutions are shown as red long-dashed lines for 
CoRot-7b, black dotted lines for Kepler 10b and blue solid lines for 55 Cnc-e. In addition, 
two solutions illustrating magnetic effects for Kepler 10-b are shown as a yellow dash-dotted 
line {Tdrag = 10"^ s) and a green short-dashed line {T^rag = 5 x 10^ s). Finally, a thick black 
solid line in Figure 2 shows the ground temperature in the Kepler 10-b model, for comparison 
with the atmospheric temperatures. 

Ignoring magnetic effects, our atmospheric solutions for the three hot super-Earths 
share many similarities. CoRot 7-b is essentially a cooler, more tenuous version of Kepler 
10b and 55 Cnc-e. This similarity stems from the strong control that vapor saturation 
equilibrium exerts on the atmospheric flow, via the exponential dependence of pressure on 
surface temperature, with a comparatively small influence for other model parameters such 
as planetary radius or surface gravity. 

In these atmospheric flows, winds are accelerated by conversion of thermal and gravi- 
tational potential energy into kinetic energy. Supersonic speeds are reached at 9s ~ 40 deg. 
Atmospheric pressures remain close to the local saturation values, within tens of % or so, 
over most of the flow, with stronger deviations in the end region at > 80 deg. Distinct 
regions of atmospheric sublimation {E > 0) and condensation {E < 0) roughly correspond 
to the subsonic and the supersonic regions of the flow, respectively. Beyond 9 ~ 85-90 deg, 
where the atmosphere becomes very tenuous, viscous drag becomes very strong and the flow 
quickly transitions into a regime where kinetic energy is returned into thermal and gravi- 
tational energy. Large Mach numbers are achieved before the atmospheric flow effectively 
comes to a halt. These flow properties are generally consistent with the ones described by 
IngersoU et al. (1985) for a sulfur dioxyde frost atmosphere on lo, despite the vastly different 
thermodynamic regime realized on hot super-Earths. Therefore, we expect these results to 
hold rather generally for hot super-Earths, including for different atmospheric compositions, 
even if a reduced angular extent of the flow (limited by strong viscous drag) can be expected 
for less massive atmospheres. 



•^Although 55 Cnc-e may require a surprisingly low-density, Moon-like bulk composition to be considered 
purely rocky (Demory et al. 2011), we include it in our analysis as a useful probe of the super-Earth 
parameter space. 
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Magnetic effects complicate this picture significantly. Magnetic drag slows down the 
winds in the subsonic region, as expected, but larger wind speeds are ultimately reached in 
the supersonic region. This delayed acceleration can be understood as resulting from the 
additional thermal and gravitational energy made available to the flow by significant ohmic 
heating, especially around and beyond the sonic point (compare the temperature profiles of 
models with and without magnetic effects in Fig. 2). The atmospheric solutions with uniform 
values of T^rag = 10^ and 5 x 10^ s for Kepler-lOb differ significantly from each other and the 
unmagnetized case, showing that magnetic effects can have a strong impact on the flow. One 
should remember that our treatment of these effects remains greatly idealized: magnetic drag 
may preferentially act on the zonal component of the flow for a dipolar planetary field and 
the magnitude of this drag (and ohmic dissipation) will vary exponentially with atmospheric 
temperature (e.g., Mcnou 2011). The large range of temperatures realized on hot super- 
Earths suggests a strongly magnetic atmospheric flow near the substellar point and a largely 
unmagnetized flow far from it (where T <^ 1000 K). These important considerations are not 
addressed by our models. 

Our solutions justify a posteriori the use of a fixed surface temperature profile, controlled 
by radiative fiuxes. The evaporation rate peaks at a few 10^^ m^^ s^^ at the substellar point, 
drops to zero near the sonic point and reaches minimal values of minus a few 10^^ m~^ s~^ 
around ^ ~ 60 deg. Given a latent heat of sublimation Cg^^b ~ 10^-10^ J kg~^ for atmospheric 
compositions of interest (e.g., Valencia et al. 2010), we estimate deviations from the radiative 
temperature profile due to latent heat effects which are at most ~ a few % at the substellar 
point. They can be neglected in view of other sources of uncertainties in the model. 

4. Discussion and Conclusion 

According to our solutions, the range of atmospheric temperatures and pressures realized 
on the day-side of a hot super-Earth is vast. Rather than a pure composition like sodium, 
it is likely that such an atmosphere would also include minor constituents which are subject 
to condensation into clouds, as discussed by Schaefer & Fegley (2009). The diverse range of 
thermodynamic conditions in the atmosphere across the dayside opens the possibility for the 
formation of many different classes of clouds. If some of these clouds are thick and reflective 
enough, they could dominate the planetary albedo properties of a hot super-Earth, given 
surface albedos of only ~ 0.1 (Legcr ct al. 2011). This may lead to interesting signatures 
in the reflected light from a hot super- Earth, if for instance strong albedo variations occur 
across the planetary dayside. A large enough albedo from clouds (> 0.1) would also generate 
variations in the surface irradiation pattern, which would imprint an additional signature in 
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the dayside thermal emission of the planet. Excitingly, the measurable optical phase curve 
of Kepler 10b could eventually reveal such signatures (Batalha et al. 2011). 

Our results show that, even for a relatively massive sodium atmosphere, the powerful 
winds which carry sublimating mass away from the substellar point on a hot super-Earth 
are unable to establish a significant atmosphere beyond the planetary day-night terminator. 
This conclusion impacts atmospheric erosion scenarios. Traditional assumptions of spherical 
symmetry and of an atmosphere at rest are poorly justified for hot super-Earths. While 
magnetized wind scenarios favor atmospheric loss from the polar regions (e.g., Adams 2011; 
Trammel et al. 2011), where there can be little atmosphere according to our models, it can 
also be the case that atmospheric temperatures in the polar regions are low enough that mag- 
netic coupling can be neglected. Furthermore, at surface pressures ^ 10~^-10~^ Pa, which 
are reached near the terminator in our solutions, the most tenuous part of the atmosphere 
may become transparent to UV irradiation (Valencia et al. 2010), which would presumably 
stall erosion. 

We conclude by noting that the sublimation/condensation rates obtained in our solu- 
tions suggest the possibility of substantial planetary resurfacing by atmospheric transport. 
The global atmospheric circulation carries mass from the substellar region to a region cen- 
tered around 9 ~ 50-70 deg, where net condensation occurs. Over Gyr timescales, this 
process could impact the atmospheric composition by preferentially transporting some (sub- 
limating) constituents relative to others and, ultimately, coupling the availability of volatiles 
at the substellar point to internal geodynamic redistribution processes. In principle, strong 
resurfacing may also influence the tidal response of the rocky planet by gradually modifying 
its gravitational moments, if condensation occurs beyond the angular extent of the dayside 
magma ocean (e.g., Leger et al. 2011). Detailed calculations will be needed to evaluate the 
magnitude of these effects. 

KM thanks Peter Goldreich for bringing to his attention the work of IngersoU et al. on 
lo. This work was supported in part by NASA grant PATM NNX11AD65G. 
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Parameter 


Kepler- 10b 


Corot-7b 


55 Cnc e 


Stellar effective temperature, (K) 


5630 


5250 


5370 


Stellar radius, (m) 


7.3 X 10^ 


5.7 X 10^ 


8.0 X 10^ 


Orbital distance, D (m) 


2.5 X 10^ 


2.6 X 10^ 


2.3 X 10^ 


Planetary radius, r (m) 


9.0 X 10^ 


10.1 X 10'5 


12.8 X 10^ 


Planetary surface gravity, g (m s"^) 


22.2 


27.1 


19.6 


Planetary substellar temperature, Tgub (K) 


3040 


2470 


2970 



Table 1: Parameters for Hot Super-Earth Models. 
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Fig. 1. — Atmospheric pressure (in Pascal) as a function of the angular distance from the 
substellar point. Blue solid, red long-dashed and black dotted lines show the results for a 
sodium atmosphere on 55 Cnc e, Corot-7b, 55 and Kepler-lOb, respectively. Solutions with 
magnetic effects for Kepler-lOb are also shown with yellow dash-dotted (Tdrag = 10^ s) and 
green short-dashed (T^rag = 5 x 10^ s) hues. Dots represent the sonic point location. 
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Fig. 2. — Temperature (in K) at the base of the atmosphere as a function of the angular 
distance from the substellar point. Same notation as in Fig. 1. A black thick solid line also 
traces the surface temperature profile adopted in the Kepler 10b models. 
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Fig. 4. — Mach number of the atmospheric flow as a function of the angle from the substellar 
point. Same notation as in Fig. 1. 



